www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第四讲 随机数的生成及随机变量抽样/liti9.m
mm=10000 Randnum1=unifrnd(-1,1,1,2*mm); Randnum2=unifrnd(-1,1,1,2*mm); xRandnum=zeros(1,mm); yRandnum=zeros(1,mm); s=Randnum1.^2+Randnum2.^2 ii=1;jj=1; while ii<mm if s(1,jj)<=1; xRandnum(1,ii)=Randnum1(1,jj); yRandnum(1,ii)=Randnum2(1,jj); ii=ii+1; end jj=jj+1; end [Y,JJ] =sort(xRandnum); plot(xRandnum(JJ),yRandnum(JJ),'.')